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AFFINE PROCESSES ON x M" AND MULTIPARAMETER TIME 

CHANGES 


MA. EMILIA CABALLERO, JOSE LUIS PEREZ GARMENDIA, AND GERONIMO URIBE BRAVO 

Abstract. We present a time change construction of affine processes with state-space 
K.™ X R". These processes were systematically studied in [DES03] since they contain 
interesting classes of processes such as Levy processes, continuous branching processes 
with immigration, and of the Ornstein-Uhlenbeck type. The construction is based on 
a (basically) continuous functional of a multidimensional Levy process which implies 
that limit theorems for Levy processes (both almost sure and in distribution) can be 
inherited to affine processes. The construction can be interpreted as a multiparameter 
time change scheme or as a (random) ordinary differential equation driven by discontin¬ 
uous functions. In particular, we propose approximation schemes for affine processes 
based on the Euler method for solving the associated discontinuous ODEs, which are 
shown to converge. 


1. Introduction 


Affine processes on the state-space E = M'” x M” are a class of processes introduced in 
[DFS03] for two reasons. First, they contain important classes of Markov processes like 
Levy processes, (multi-type) continuous branching processes with immigration, and of 
the Ornstein-Uhlenbeck type. That is, they contain the fundamental examples of models 
in (stochastic) population dynamics (as in [Lam08]) and mathematical finance (as has 
been argued in [DFS03] and [Kal06]). Second, they are analytically tractable. Indeed, 
they have been shown to be parametrized in a manner similar to Levy processes and 
one can access their finite dimensional distributions by solving an ordinary differential 
equation of the Riccati type (cf. [DFS03]). 

To define them, let Z = (Z,,f > 0) denote a stochastic process on a measurable space 
(Q,^) whose paths are cadlag functions from [0,°°) to E. The canonical filtration of 
Z will be denoted Suppose that the measurable space is equipped with a family of 
(sub)probability measures (IP^jZ E E) such that under each the process Z starts at z. 
Furthermore, we assume that Z is stochastically continuous under P^ for any z 6 E and 
that these measures constitute a Markov family: 


E,(/(Z,+,) I ,F;) = Prf(Zs) where PJ(z) = E,(/(4)). 


Definition. The Markov family (P^jZ E E) is affine if 

( 1 ) 



for all u E x /M”, where M'” denotes the set of elements of whose coordinates 
are negative. 
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These proeesses are part of a larger one of so-ealled affine proeesses on general state- 
spaees and mueh reeent work has been aimed at eharaeterizing these Markov proeesses, 
for example by proving their Feller property and the preeise form of their infinitesimal 
generator. This work started in [DFS03] for regular affine proeesses on E and was 
later extended in [KRSTll], [KRST13] and [CT13] by proving that regularity already 
follows from stoehastie eontinuity and also by eonsidering more general state-spaees 
than we do here. 

Our main result aims at giving a pathwise eonstruction of affine processes in terms of 
a multiparameter time change of Levy processes, which are considered as more basic 
building blocks. 

Theorem 1. LetX ^,... and Y be independent Levy processes on M™+". We suppose 
that the first m coordinates ofY are subordinators, that X‘^' has no negative jumps and 
thatX^'J is a subordinator for 1 <i,J <m and i j. Furthermore, in the Gaussian part 
ofX\ the i-th coordinate is assumed independent of coordinates m + 1 up to n. 

Let fi be an n X n real matrix. Then, for any z 6 E there exists a unique solution Z to 



I < j <m 
Pij i<j<n 


( 2 ) 


with 



whose first m coordinates are non-negative. IfF^ denotes the law ofZ, then 6 E) 
is an affine Markov family on E and every affine Markov family E is obtained by this 
construction. 

Note that the non-negative coordinates are more difficult to handle. Indeed, the non¬ 
negative coordinates alone constitute an affine process with n = 0 , which is then called a 
multitype continuous-state branching process with immigration (CBI) introduced (with¬ 
out immigration and with m = 2) in [Wat69]. Once we analyze the case n = 0, we will 
then get the general case by solving a linear differential equation driven by the solution 
when n = 0. These real-valued coordinates constitute the Ornstein-Uhlenbeck part of 
the process, which is now not only driven by a Levy process but also by a sum of time 
changed Levy processes. 

Equation (2) represents a multiparameter time change equation proposed in [KurSO] 
to generalize the classical time change construction of Markov processes of Volkonskii 
(cf. [Vol58], [Dyn65, Vol 1, Ch 10]). A multiparameter time change representation 
of affine processes was first proposed (in a weak sense) in [Kal06]; in that paper, the 
question of whether the affine process was adapted to the filtration of the Levy pro¬ 
cess was left open. Recently, there have been a number of results concerning this time 
change representation. For example, the PhD thesis [Gab 14] (the relevant chapter is 
found in [GT14]) proves existence (under additional but minor technical assumptions) 
for a time change representation as in Theorem 1 . A discrete space version of Theorem 
1 has also been recently studied. Indeed, a construction of Galton-Watson processes 
(without immigration) in terms of multiparameter time changes of random walks is 
found in [CL13] in discrete time and [Chal4] in continuous time. More generally, the 
connection between time changes and changes of measure and the application to math¬ 
ematical finance is explored in [BNSIO]. The main contribution of our work is that we 
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prove uniqueness of the pathwise representation in (2) (as well as for an accompany¬ 
ing inequality). Uniqueness is the main tool in the forthcoming stability analysis of the 
pathwise representation. 

We now state some continuity properties of the system of equations (2). Consider 
a sequence, indexed by / > 1, of m-|- 1 stochastic processes which 

satisfy the upcoming hypothesis H of p. 6. Consider also any sequence of numbers 
0 < (7/ ^ 0. The number C/ is interpreted as the discretization parameter to be used 
in an Euler type scheme as follows. When Oi > 0, let and C■'’^ 1 < j < m + n,he 
defined recursively by means of 


(3) cy=0, 


zi’^, = 






i=\ 




when 1 < 7 < m, while for 1 < y < n we only change the definition of to 


(4) 


7m+j,l 

"CTik 


= z/ + £r’'"+^4oc; 

i=l 


i,l ,ym+j,I 
Oik ' ^ Ojk 


+ Ec: 

1=1 


m+i,l 
Oik 


A. 


When (7/ = 0, the forthcoming Lemma 3 asserts that (2), when driven by 
admits a (global) solution (which could, in principle, explode). In that case, we let 
be any such solution. We recall in Subsection 6.2 the definitions of the Skoro- 
hod Ji topology and of the uniform Ji topology. 


Theorem 2. LetX^ ^and Y be as in Theorem 1. LetZ^C be the unique processes 
satisfying (2). 

Suppose thatX^^'^\ ... are stochastic processes which satisfy hypothesis H 

of p. 6 and such that X^’ ’^ converges to X^’' (and 7 ’^ converges to Y) as I ^ (The 
convergence can be weak or almost surely in the Skorohod J\ topology when (2) has no 
explosion and in the uniform Ji topology in case of explosion. ) Assume that zj —)■ zT 

If Z '^, C’^ are any processes satisfying (3) and (4) when Oi > 0 or (2) with respect 
to the driving processes . ,X'^’'’^,Y'’^ when Oi = 0 then —>■ C (with respect 

to the topology of uniform convergence on compact sets when there is no explosion 
and pointwise in case of explosion) and Z'’^ —)■ Z' for I < i < m (with respect to the 
Skorohod Ji topology when there is no explosion and in the uniform J\ topology in case 
of explosion) as I ^ (The convergence will be either weak or strong depending on 
the type of convergence of the and Y.) 


Note that the above limit theorem is either weak or strong, which follows from con¬ 
tinuity properties of the multiparameter time change equations explored in Section 6. 
Indeed, we believe this is one strength of the time change representation versus, for 
example, the SDE representation which is found in the one-dimensional case in [ELIO] 
and [Li 14]. Indeed, even in the case of continuous sample paths, it is known that solving 
SDEs is a discontinuous operation of the driving processes. A manifestation of this is 
found in Wong-Zakai type phenomena (discovered in [WZ65]) and depending on the 
type of approximation to the driving processes one obtains limits to different SDEs, as 
has been argued in [EH 14]. On the other hand. Theorem 2 does not depend on how 
one approximates the driving processes. We are not advocating, though, the use of one 
representation over another. The construction of [Li 14] is useful in the genealogical in¬ 
terpretation of continuous branching processes, constructing directly some of the flows 
in [BLG03]. 
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From Theorem 2 we deduee a limit theorem eoneeming multi-type Galton-Watson 
proeesses stated as Corollary 1. In the one-dimensional case, Corollary 1 includes limit 
theorems found in [Cri74], [Li06] and [CPCUB13]. The multidimensional case has of¬ 
ten been studied in the literature when the limit process is continuous, as in [JM86]. 
We state a version without immigration, just to illustrate the kind of statement one 
can achieve as well as the technique. The technique can be adapted to the case of 
immigration as in Corollary 7 of [CPCUB13]. Let 1 < / < m) be independent 

J-dimensional random walks. Suppose that has jumps in Z greater than —1 and 
that otherwise the coordinates have jumps in M. Let = (k|,..., G N™ be a se¬ 
quence of starting states and define recursively the sequences I < j <m) 

and 1 < y < m) by 


^0 


= 0, Z/, = k‘ 


i=l 


and C 


7 

n+l 


= cl+zi 


'n+l' 


It is easy to see that for each Z, Z^ is a multitype Calton-Watson process such that the 
quantity of descendants of type j of an individual of type i has the same law as 
when i ^ j and the law of -|- 1 in the remaining case. However, if X^ is extended 
by constancy on intervals of the form [n,n-\-\) with n G N, we see that is the Euler 
type approximation of span 1 applied to X^ that we have just introduced and Z^ is the 
right-hand derivative of 


Corollary 1. Let X^’ ’^, I < i < m be independent d-dimensional random walks. Sup¬ 
pose that has jumps in Z greater than —1 that otherwise the coordinates have 
Jumps in N. 

Assume that for each i in {1,... ,ni} there are scaling constants ai and b\ for Z > 1 
such that 

—Xy’V>0,l < j<m 

bj b‘t ^ ^ -■> - 

converges in Skorohod space (either almost surely or in distribution) to a Levy process 
X'’'. Furthermore, a/ —)■ oo, Z?^/a/ —)■ oo and kj is such that kjai/bj -+ zL 
Then, the scaled Galton-Watson processes 

>0,1 <j<m 

started from 1 < 7 < m) converge in Skorohod space (either almost surely or in 
distribution) to the unique CB process Z started from z and constructed from X and 
Y = 0 in Theorem 1. 


We end this section with an application of Corollary 1 . Note that the different pro¬ 
cesses in Corollary 1 have scalings that have to be adequately balanced in order to obtain 
a limit (with non-trivial reproduction and immigration components). In order to exem¬ 
plify how this could be done, let us start by considering the framework of Theorem 4.2.2 
of [JM86], giving a limit theorem for nearly critical multitype Calton-Watson processes 
under finite-variance assumptions. Indeed, consider a sequence of multitype Calton- 
Watson processes Z ’^ such that is the law of the offspring of an individual of type Z. 
We then define the mean matrix M by means of 

<■= E 
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Assume that = Id +C //1 where Ci ^ C as I ^ oo. Consider also the varianee matrix 

given by 






1/2 


Supose that ojj —)■ 0 ';l,=y as Z —)■ 0 ° and that the following Lindeberg condition holds: 


ki>e^/n 


as / —)■ 00 . Recall our construction of such a process in terms of random walks 
for i = 1,..., m. From our hypotheses, it follows that X'^ 2 ’^/^ converges to OiB' + Ci^i Id 
where 5' is a standard Brownian motion. Indeed, the convergence of one-dimensional 
distributions is deduced from the Lindeberg-Feller central limit theorem. Because of 
independence and stationarity of the increments this implies the convergence of finite¬ 
dimensional distributions and tightness is easily deduced from the Aldous criterion. For 
i 7 ^ j, one sees that /1 Cij Id as / —)■ 0 °. Indeed, it suffices again to establish con¬ 
vergence of one-dimensional distributions which follow from Chebyshev’s inequality. 
Tightness again follows from the Aldous criterion. Hence, Corollary 1 allows us to con¬ 
clude that if Zq/Z —>■ z then Zj /I converges weakly to a continuous branching process Z 
with continuous sample paths. One can then use the martingales associated to X, as in 
[KurSO], to see that the generator of Z is given by 


m 


E 


Ziof 

2 dzj 





This fact can also be deduced from the infinitesimal parameters of Z that are introduced 
in Section 2 and from the proof of Theorem 1 . 

Our work continues and extends the one-dimensional situation covered in [CPGUB 13]. 
There are however, important differences with that work. First of all, the discussion of 
uniqueness to (2) now relies on the concept of (lack of) spontaneous generation. This 
is to be contrasted to the previous analysis based on taking inverses. The multiple time 
changes make this one-dimensional approach unfeasible. On the other hand, we also 
take the point of view of multiparameter time changes from [KurSO], providing a very 
concrete (but general) example of its applicability. This has led to several simplifications 
when proving that solutions to ( 2 ) are affine processes. 

The paper is organized as follows. We first consider a deterministic framework for 
equation ( 2 ) when n = 0 and analyze existence, uniqueness, and basic measurability 
questions. This is done in Section 3. We then undertake the proof of Theorem 1 when 
n = 0, which reduces basically to establishing the Markov property and constructing 
relevant martingales, in Section 4. The case of general n is taken up in Section 5. 
Finally, we pass to the stability of equation 2, which contains the proofs of Theorem 2 
and Corollary 1 in Section 6 . 


2. Preliminaries on affine processes 

Let Z be an affine process with laws G E). Let d> and i// be defined as in Equa¬ 
tion (1); applying the Markov property, we get the semi-flow property 

(5) ^/(t-l-ijM) = V/(5, i//(Lm)) and ^{t + s,u) . 
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From Theorem 5.1 in [KRSTll] or Theorem 3.3 in [KRST13], it is known that the 
following derivatives exist and are continuous as a function of u: 




f=0 


and R(u) = —\j/(t,u) 
ot 


r=0 


From the semi-flow property, we deduce the so called Riccati equations 

( 6 ) 

with the initial conditions 


—and —\j/{t,u)=R(\j/{t,u)) 
ot ot 


^>(0,m) = 1 and V^(0, w) = u. 

If were non-zero and (j) were continuous and satisfied e‘^ = 4>, we would obtain the 
more familiar equation 

which gives 

^{t,u)= / F{\{f{s,u)) ds. 

Jo 


Furthermore, Theorem 2.7 in [DFS03] asserts that F and R have the following very 
specific form: if X \ . • •; and Y are Levy processes satisfying the conditions of The¬ 
orem 1 then F and R= (Ri,... ,Rm+n) are the unique continuous functions such that 

(7) and W = E ) 

for 1 < / < m while for m -|- 1 < / < n we set 

n 

Rm+iitJ) — 

7=1 


Furthermore, Section 6 of [DFS03] discusses the (global) existence and uniqueness 
of the generalized Riccati equations of (6). 


3. PATHWISE ANALYSIS OF THE MULTIDIMENSIONAL TIME CHANGE EQUATION 

Following [CPGUB13], we begin by considering a deterministic system of time 
change equations appearing in Theorem 1 in the case of non-negative processes (n = 0). 
Consider m{m+\) cadlag functions labeled {f‘J, 1 < i,j < m} and {g-', 1 < y < m}. 
These functions satisfy the following requirements: 

HI: f’J has no negative jumps if i = j and is non-decreasing otherwise. 

H2: is non-decreasing. 

H3: g^(0) -f 1 fJ{0) > 0 for 1 < 7 < m. 

The above hypotheses are collectively denoted H. 

We seek a solution to the following system of equations for the cadlag function h = 

(/ii, 

.ft . 

(8) h-^{t) = oc\t)+g-^{t) for 1 < 7 < m where c-’{t)= / h\s)ds. 

Jo 

This system can also be thought of as an ordinary differential equation for c when one 
notes that is the right-hand derivative of cL With this interpretation, we might want 
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to use other initial eonditions for c rather than only zero. This amounts to shifting the 
funetions note however, that the shifts must still satisfy H3. 


3.1. A basic monotonicity lemma and existence. Our approach to the study of (8) is 
based on its monotonicity properties. We begin with a simple and useful case of this 
and postpone an elaboration of this idea which will be useful to obtain uniqueness. 

Lemma 1. Suppose that we have two sets of functions P = and P = 

satisfying hypothesis H. Assume that p’J < p’J and g-^ < g^ and that additionally, for 
every y G {1,..., m}, either or gJ < gh Ifh and h are non-negative functions 

that satisfy (8) driven by P and P respectively then c <c. 


Proof Let 

P = (/’■' ,g^A< i, j < m) and P = {f'fgf 1 < i, j < m) 


be a system of functions satisfying the assumptions of our lemma and h, h the associated 
non-negative solutions to (8). For any e > 0 and any a > 1, define c^{t) = c^{e -\- at). 
Hence W has a cadlag right-hand derivative h'’ given by h'’{t) = ahpe -f at). We then 
define 

T = inf {t > 0 : c-’f) > c^{t) for some _/} 

as well as the set J of indices y G {1,..., m} such that c^ exceeds W strictly at some 
point of any right neighbourhood of T. If y G 7 then c^ (t) = C“' (t) while c'(t) < c'(t) 
for i f j, so that also oc'(t) < P’j o c'(t) for i j. We deduce the following for 
jeJ: 


0 < = ^/'’■' oc'(t) +g-’{T) 

i 


i 


< a 


^/'’■'oc'(T)+|^(e + aT) 



(Note that the right-hand side of the first strict inequality cannot be zero, which justifies 
the second strict inequality.) We deduce that c-’ remains below c-' in a right neighbour¬ 
hood of t which contradicts the definitions of T and J. We deduce that c^ < c^ and, 
letting a go to 1, that c^ <cf □ 


We now tackle existence for (8) in the case when only p’f I < j <m are not piece- 
wise constant. The proof will be based on the observation that under the piecewise 
constant hypotheses, the system (8) is one-dimensional on adequate intervals. The 
piecewise constant case will allow us to prove existence for (8) in general through the 
monotonicity proved in Lemma 1 . 

Lemma 2. Let {/'’Z gf 1 < f, y < m} satisfy H and suppose that and gj are piece- 
wise constant i/l < /, j<m and i f j. Then, there exists a solution h = \h^ : 1 < y < m} 
to (8). This solution exists on an interval [0, t) and c\_ = oofor some j. 

The time T is termed the explosion time of c. 
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Remark 1. For the one-dimensional case, existence follows from Theorem 1 in [CPGUB1 3] 
which asserts that the problem IVP(/,0,Jc) consisting of a finding a function c with a 
right-hand derivative h which satisfies 

IVP(/,0,Jc) : h = foc with c(0)=t 


admits, for any x >0 and any cadlag function / such that / has no negative jumps, a 
unique solution which lacks spontaneous generation. When f{x) = 0, the only solution 
lacking spontaneous generation is the function c{t) = x. When f{x) > 0, the unique 
solution can be constructed by a Lamperti type transformation obtained by first making 
zero absorbing after x; formally 


r = inf{t > X :/(t) = 0} and /(t) = 
We then define i on [x, oo) by means of 




t < T 

re [r,oo] • 


Note that i is strictly increasing on [x,T) and infinite on (T,oo). Then, let c be the right- 
continuous inverse of i (in the sense of Lemma 0.4.8 of [RY99]). Note that c is strictly 
increasing on [0,/(T—)] and constant on [/(T—),oo] and by definition c(0) = x. Then, 
since the right-hand derivative of i exists and equals 1 //, then c also admits a right-hand 
derivative (on [0,/(T—))), say/?, and we have/i= l/(l//oz^^) =foc. The function c 
so constructed from / is called the Lamperti transform of / absorbed at its first zero after 
X . Note that c(oo) = T. In the one-dimensional setting, whenX is a spectrally positive 
Levy process. Proposition 2 of [CPGUB 13] shows that there is a unique solution C to 
IVP(x-|-Y,0,0) (with right-hand derivative Z) which has zero as an absorbing state; if 
T denotes the hitting time of zero of x + X, X equals X stopped at T, then C is also 
the unique solution IVP(x-|-X,0,0), so that Coo = T. This one dimensional result is 
important in our proof of uniqueness of solutions to 2. Since stopping a cadlag process 
at a stopping time and looking at a cadlag process at a random time are measurable 
transformations, we see that the Lamperti transformation is measurable on the Skorohod 
space of cadlag trajectories with the c-field generated by projections. This would hold 
even if we take the initial value x to be random and measurable. 


Proof. Suppose first that 


f‘-‘ = E 




k=\ 




= 


HyiK, 


k=\ 




where x'f’_^ < x''^’ if i 7 ^ 7 , 0 = Iq-’ < < ■ • •, the sequence k>0 has no accumula¬ 

tion points (similar assumptions hold for gj) and additionally, for each j 

f’J{0) + Y,xY+y{>0 

¥j 

so that assumptions H hold. Let (resp. Tj) denote the set of change points of the 
functions (resp. g^): 

r-J = {'f : 0 < tj and r;={<L0<*}. 
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Let To = 0 and, for any 7 = 1,..., m, let Cj be the unique solution of the problem 
IVP , 0,0 j, where the function f( is given by 

¥j 

We now define the times 

tJ = inf {t > 0 : t G 7]- or c\ (t) G } and Ti = min tJ . 

Set c-i equal to c{ on [0, Ti] and recursively define ^ as the solution to IVP j , 0, c-^ (T„) 

where the function is given by 

fn+li^) = °c'(Tn) +g\Tn) ■ 

¥.i 


We then define 


= inf {t > T„ : t G or c\ {t - T„) G } and T„+i = niin T^_^i 

and let ¥{t) = — t„) on [t„,t„+i]. We assert that c = (cV--:C™) solves ( 8 ); 

the proof is by induction. However, note that the starting point of c„ is chosen so that 
c is continuous and has a cadlag right-hand derivative. On [0, Ti], o d and g^ are 
constant and hence, equal to their value at zero. Hence, if we let y stand for the right- 
hand derivative of d, we obtain the following equalities for any t < Ti 

*■'(<) =// 

¥j 

¥j 

which allow us to conclude that c solves ( 8 ) on [0, Ti]. On the other hand, if we assume 
that c solves ( 8 ) on [0, T„], then note that, by definition, o d and g-^ are constant on 
[t„, t„+i] . We deduce that for t G [t„, t„+i] : 

= f’-’ °c^+i ° +g-'o c ^( t „) 

¥j 

=° c' (t) -F (t) 

¥j 

so that c solves ( 8 ) on [ 0 , T„+i]. 

Since T„ increases in n, there are two possibilities: either T„ —)■ 0 ° (in which case the 
solution we have constructed is a global solution) or c-^ (t„) —)■ 0 ° for some j by definition 
of Tn, T„, and the fact that the sets Tij and Tj have no accumulation points. In the latter 
case, c explodes. □ 


Remark 2. As in Remark 1, we note that if we apply the procedure of the above proof 
in the case of cadlag stochastic processes (satisfying the conditions of Lemma 2) then 
the solutions are measurable. This follows because on adequate intervals (which are 
obtained by stopping), the solutions are unidimensional and are constructed through the 
Lamperti transformation. 
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We now tackle existence for (8). 

Lemma 3. Let 1 < < m} satisfy H. Then, there exists T > 0 such that a 

non-negative solution h to (8) exists on [0, t). Furthermore this solution explodes at x 
and is maximal: 

(1) If c^{t) = f^y^s) ds then c-^(t—) = oo for some j. 

(2) Ifh is another solution to (8) (with its corresponding c) then c <c on the interval 
of existence ofh. 

Proof For I <ij <m with i j consider a sequence of cadlag functions fy and gi 
which are piecewise constant, are strictly bigger than and gf and decrease as n —)■ oo 
towards and gl respectively. We then set ff^ = fif Using Lemma 2, we can 
consider for any n a solution = {h\, ..., h™) to (8) driven by {fh\gi} . By Lemma 1, 
we see that the cumulative population of exceeds the cumulative population of any 
solution to (8). 

Fix any K> 0 and use it to stop c„ at the instant Xn^K that any one of its coordinates 
reach K. Call the resulting function c„. Since c„+i < Cn then Xn,K < T:n+i,K', set Xk = 
lim„ Xn,K- Note that Cn has a cadlag derivative h" given by 

Kit) = h<T„^KKit) ■ 

Hence, hi can be bounded on any interval [0,t] by mmaxijmfx<K fK(x) +giit), and 
can then be bounded in n by construction of fK and gi. By the Arzela-Ascoli theorem 
(which applies since c/,(0) = 0), c„ is sequentially compact. We now show that every 
subsequential limit coincides. Indeed, if c is the (uniform) limit (on compact sets) of ci/^ 
as k —)■ oo, then the bounded convergence theorem implies that for any t < Xk- 

y{t)=\imcl{t)=lim f °Kis) + ds = [ °^'i^)+Kis) ds. 

We conclude that c admits a right-hand derivative h on [0, Xk) which satisfies (8) on 
[0, Xk) . However, c is the maximal solution by construction (since we can apply Lemma 
1 to the approximations c„), so that all subsequential limits agree on [0,%]. Finally, 
note that before Xk the coordinates of c have to be smaller than K and that at Xk some 
coordinate equals K. Hence Xk coincides with the instant in which some coordinate of 
c reaches K. By uniqueness, one can construct a function c which coincides with c on 
[0, Xk), so that c is defined and solves (8) on [0, x) where x = limA- Xk- By construction, 
c explodes at x and is maximal in the class of solutions to (8). □ 

Remark 3. Recall that the approximations of the above proof are measurable in the case 
of applying them to cadlag stochastic processes thanks to Remark 2. Then, applying the 
construction to a cadlag stochastic process X satisfying hypotheses H, we get another 
pair of stochastic processes Z and C. Since Z and C are cadlag, then X‘d oC‘ is also a 
stochastic process. 

3.2. Spontaneous generation and minimal solutions. An interpretation for the one¬ 
dimensional case of (8) was proposed in [CPGUB13] by noting that if represents 
the breadth-first walk on a (combinatorial) forest representing the genealogy of a pop¬ 
ulation with immigrants along each generation and codes the immigration to the 
population then K is the population profile (that is, the sequence of generation sizes), 
while is the cumulative population. The multidimensional case of this discrete cod¬ 
ing can be found in Subsection 2.2 of [CL13], when gl = 0 for all j, and it shows that 
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the one-dimensional interpretation still holds. In partieular, the diserete interpretation 
gives sense to the following definition of laek of spontaneous generation: in the one 
dimensional ease, a solution h = f o c + g laeks spontaneous generation if h{s) =0 (the 
population is zero at time s) and g is eonstant on [ 5 , 1 ] (there is no immigration on [s,t]) 
implies that h = 0 on [ 5 , 1 ] (the population remains at zero). Perhaps it is then surprising 
that there are solutions featuring spontaneous generation, but when / is the typieal path 
of a normalized Brownian exeursion then they do exist (cf. [CPGUB13, Sect. 2]). 

Definition. Let ,g-^) satisfy H. We say that a solution h = to (8) has no spon¬ 
taneous generation if whenever ( 5 ) = 0 for some 5 > 0 and for all 7 in 7 C {1 ,..., m} , 
we have that the strict increase of at s for some j E J implies that either gj increases 
strictly at s or there exists i ^ J such that o c‘ increases strictly to the right of s. 

As a remark, we mention that h lacks spontaneous generation if and only if at any 
5 > 0 such that the set J{s) = {j '■ h^{s) = O} is nonempty, the strict increase of at 5 
for some j G J{s) implies that either gj increases strictly to the right of s or there exists 
i ^ J{s) such that increases strictly to the right of c'(5). 

The definition works very well with induction on the dimension, in the sense that 
if h = {h\i < m) is a non-negative solution to (8) driven by (/'’^, 1 < i,j < m) and 
{g-’ ,j < without spontaneous generation and mi < m, we can then consider 
as a solution, which will lack spontaneous generation, to (8) but driven by and 
g-’ + Zi>mi ° c‘ for 1 < 7 < mi. 

The importance of solutions lacking spontaneous generation is that they have mono¬ 
tonicity properties (see Lemma 4 below) and, consequently, they are minimal solutions 
to (8) as well as unique. In particular, if all solutions of (8) can be shown to have no 
spontaneous generation, then there is at most one solution. There are two cases when 
we can actually apply this technique. First, when g-’ is strictly increasing for all j since 
then solutions trivially have no spontaneous generation. Another example is when (8) 
is driven by Levy processes satisfying the hypotheses of Theorem 1 : we will show in 
Lemma 6 that solutions have no spontaneous generation, which covers the uniqueness 
statement in Theorem 1 . 

Lemma 4. Suppose that we have two sets of functions P = {f'fg^) and P = {f'fg^) 
satisfying hypothesis H. Assume that f’j < and gj < gf Ifh and h are non-negative 
functions that satisfy (8) driven by P and P respectively and h lacks spontaneous gen¬ 
eration, then c < c. Hence, (8) admits at most one solution h whose coordinates are 
non-negative and have no spontaneous generation. 

The above lemma also tells us that solutions without spontaneous generation are min¬ 
imal in the sense that their primitive is a lower bound for the primitive of any other 
solution. 

Proof This proof is an elaboration of the proof of Lemma 1. We proceed by induction. 
Let m = 1, let F = (/,g) and P = (/,g) satisfy hypothesis H and let h and h be non¬ 
negative solutions to (8) driven by P and P and lacking spontaneous generation. Let 
£ > 0 and a > 1 and use them to define c by means of c{t) = c(£ -f at). Let 


T = inf{t > 0 : c{t) > c(t)} . 
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Note that c has a cadlag right-hand derivative h given by h{t) = ah{£ -h at). If /i(t) > 
0, sinee c(t) = c(t), our assumptions give 

h{T) =foc{T)+g{T) 

</oc(T)+g(T) 

< a [/oc(t) -i-|(e + aT)] =h{x). 

Henee, c < c on a right neighbourhood of T eontradieting its definition. On the other 
hand, if h{x) =0, then we ean only infer, as in the previous display, that 

0<h{x) = foc{x)+g{x) 

<foc{x)+g{£ + ax) 

< foc{x)+g{£ + ax) = 0. 


We eonelude that g is eonstant on [t, £ +ax] whieh, by laek of spontaneous generation, 
shows that = 0 on the same interval so that c eannot exeeed c in any small enough 
right neighbourhood of x. We eonelude that c{t) < c(e -f ax) for any t >0, any e > 0 
and any a > 1. Henee c < c. 

Let m> 2. Suppose now that the monotonieity statement c < c is true for any so¬ 
lution to (8) of dimension strietly less than m. Let P = ,g-^, 1 < < m) and 

P= <hi< m) be a system of functions satisfying the assumptions of our 

lemma in dimension m and h and h the associated non-negative solutions to (8) without 
spontaneous generation. We proceed as in the one dimensional case: for any e > 0 and 
any a > 1, define C“'(t) = c-^(e -f at). Hence c-^ has a cadlag right-hand derivative h’ 
given by h’{t) = ah^{£ -f at). We then define 

T = inf {t > 0 : c\t) > c\t) for some j} 

as well as the set J of indices y G {1,..., m} such that d exceeds strictly at some 
point of any right neighbourhood of x. If j G 7 then c^{x) = C“'(t) while c'(t) < c'(t) 
for i 7 ^ j. \fh'{x) >0 for some j G J, we infer that 

h\x) = fj’-' oc-’{x) + ^f’-^oc\x)+g-’{x) 

¥j 


< a 


p,J o ¥1^ 


'+E/' 

¥j 




T)-M''(e-f oct) 


= hHx) 


We deduce that c-' remains below ¥ in a right neighbourhood of x which contradicts the 
definitions of x and J. Hence, we can assume that h\x) =0 for every j G J. Note that 
if 7 = {1,..., m} then 

0 < h-^{x) = +g-'(T) 

i 

<E/'’''°c'(T)+g-'(e + aT) 

i 

<E/'’-'oc'(T)+|-'(e + aT) =0. 


We conclude not only that hj{x) = 0, but also that gj is constant on [t, e -f ccr], which 
implies, by lack of spontaneous generation, that is constant on [t, e-f ax], which 
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contradicts the definition of T. Henee, we ean assume that 7 C {1,..., m} and by rela¬ 
belling, we write J = {1,... ,mi} where m\ < m. For every j > mi, d < W in a right 
neighbourhood of T. Also, note that , j < mi solves system (8) in dimension mi when 
driven by [f'^,i,j<mi) andg-' -t-£,>„j/'’-^ oc“'. The same remark holds for c-^, 7 < mi. 
Sinee this system have dimension strietly less than m and the redueed system hi,. 
has no spontaneous generation, monotonieity holds for them and we ean eonelude that 
c-’ < C“', j < mi in a right neighbourhood of T, whieh again eontradiets the definition of 
T. As before, we eonelude that c < c. 

Finally, suppose that h and h are two non-negative solutions to ( 8 ) whieh laek spon¬ 
taneous generation and are driven by the same funetions P. Applying our monotonieity 
statement, we see that c — c whieh then implies h = h. □ 


3.3. Further consequences in the stochastic setting. We now show that the proeess 
C is a multiparameter random time ehange in the sense of [EK 86 , Ch. 6 ]. For this, 
eonsider the a-field 


(9) 






^1 




im 




Lemma 5. Let X'” and Y be a stochastic process satisfying hypotheses H. Let Z 

be the solution to (2) (with n = 0) such that its primitive C is maximal. Then 


a{Zi,Ci,YJ,T 


Reeall that the solution construeted in Lemma 3 has a maximal primitive. The proof 
will be based on a Galmarino type test in the multiparameter setting. (Cf. [Gal63] and 
[RY99,Ex. 1.4.21].) 


Proof. Let Z, C be as in the statement. Consider also the solution Z to (2) (with m = 0), 
but now driven by X' stopped at (denoted X') for 1 < i < m, by 7 stopped at t (denoted 
Y), and sueh that its primitive C is maximal. Analysing the eonstruetion of C, Z, C and 
Z we note that if C] < ti for \ <i<m then Z = Z, X'-'^ o C' = X'-'^ o O and C = C on [0, t]. 
Sinee Z,, o O and O are measurable funetions of X^'f the statement follows. □ 


We now study the uniqueness of (2) when n = 0. 

Lemma 6. LetX^,... ,X'” and Y be Levy processes satisfying the assumptions of Theo¬ 
rem 1 when n = 0. Then, almost surely, solutions Z to (2) have no spontaneous genera¬ 
tion. 


Proof. Let Z (equivalently C) be a solution to (2). Let x be the first instant such that 
Z admits spontaneous generation. We argue that t = oo by contradietion. Indeed, we 
first show that T > 0 almost surely and then we apply arguments related to the strong 
Markov property to deduce that T < T on the set T < oo. 

Let us now show that T > 0. Note that T > 0 means that if 7 = {] <m: Zj = 0^ 
then the assumption that 7-^ = 0 on a right neighbourhood of 0 for every y G 7 and that 
A'’-' = 0 on a right neighbourhood of zero for every 7 G 7 and i ^ 7 implies that Z^ = 0 for 
every 7 G 7 on a right neighbourhood of zero. Hence, the problem is reduced to proving 
that if Zj = 0 and 7-^ =0 for every 7 then Z = 0. (This corresponds to analyzing the case 
of multitype CB proeesses without immigration. ) Let T > 0 be sueh that Cj <°° for 
every 7 . Sinee is a subordinator for any i 7 ^ 7 then /h exists and equals 

the drift eoeffieient of A (cf. [Ber96, Prop 8 , Ch III]). Hence, there exists M > 0 sueh 
that A' ’j < Mid on [0,Cy] for all i 7 ^ 7 . If 7 is any coordinate such that X-^’j is a finite 
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variation Levy process (that is, the difference of two subordinators), the same argument 
implies that (for a possibly different M) X-^’j < Mid on [0,Cy]. For this coordinate we 
see that 

0<zi <MCi +mY^C\. 

¥j 

From Gronwall’s inequality, we see that Z equals zero until there exists i 7 ^ j such that 
O grows. Hence, coordinates j such thatX-^’-^ have finite variation cannot be responsible 
for spontaneous generation. To analyze coordinates with infinite variation, recall from 
[Rog 68 ] that if has infinite variation then liminf/j^oZ^’-^/fi = — 00 . We consider 
m > 2 since the case m = 1 has been handled in [CPGUB 13]. In particular, for A > m — 1, 
we can choose a sequence (g) decreasing to zero and such that Xf^^ < —MAtn. We then 
choose e„ in the interval (Mt„,MAG/(m—1)). Now, define = £„ V (Mid) for/ 7 ^ j 
and consider a solution Z to 


z/ = x^'^ O c/ + O c\. 

¥j 

A modification of the proof of Lemma 1 shows that < & on [0, T]. However, note 
that while every & is below e„/M, & behaves as the solution to 

z/=A^'’^'oC/ + (m-l)e„. 

It follows that if Z-' reaches zero before any O exceeds e„/M then Z^ remains at zero 
afterwards (since this happens for the one-dimensional problem defining Z^). However, 
recall from Remark 1, that in the one-dimensional case the total population {CL) equals 
the time the reproduction function reaches zero (inf{t > 0 : (m — 1 )e„ = 0}). Since 

(m — 1) £„ -\-Xf^^ < (m — \ )£n—MAtn < 0 and < £n/M, it follows that Z^ reaches zero 
before & reaches Sn/M. Hence, & < Sn/M and so C = 0 on [0, T]. We have hence 
shown that T > 0. 

We now the following identity in law: 

(10) 

(x‘(c; + -)-x>(c;),...,x"(c" + -)-x"(cn.i',+,-n)|T<“>=(x‘,...,x'".y). 

The identity in law (10) implies a contradiction since if Z satisfies (8) but driven by 
the left-hand side of (10) with initial value Z^ then Z should, by definition of T, have 
spontaneous generation at time 0, which is impossible. 

To prove the identity in law of (10), we first prove that for any 0 < we have 

(11) A = {cj < 1 ,....,C? < T> 1 } e 

Indeed, the above will follow from proving that 

{T<t}e a{ci,Zi,Yj,r’JoCl :0<s<t,l<ij<m). 

thanks to Lemma 5. However, note that T < t is and only if there exists J C [m] = 

{1,..., m} and j G 7 such that Z^ presents spontaneous generation over a common in¬ 
terval of constancy of and X‘’j o C' of length greater than £. This can be discretized 
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as follows: 

p<<}= u uu u n u nn 

yc[m] i^J e>0 qe[0,t] 5>0p<q-£i^J feJ 
eeQ 5eQ peQ 

|z^ > 0,0 < ^] ’ < 51 n o C\ 7“' are constant on [p, q ]} , 

where = inf |zi : ^ e 

We now consider the random times and Q where 

Xn = {k+ l)/T ifk/T <x<{k+ \)/T 

and 

c = {h + l)/2'’ if < c;, < [h + l)/2^ 

Then, thanks to (11) 

{T=(*+l)/ 2 ",c; = (*i+l)/ 2 "} 

= [k/i- < T < (*+ i)ir,k,ir < < (fc+1)/2"} 

a 

^(fci + I)/2«,...,(fc„+I)/2«,(fc+I)/2« 

Also, note that T„ and Q decrease to T and Q respectively. 

Consider now the processes Z' and Y where X} = Z' — and 

C„4-r 

We assert now that the joint law of A \ ... ,X'" and Y equals the law of X \ and 

y. To prove this, we focus on the one-dimensional distributions since the computation 
of the finite-dimensional distributions is just notationally more cumbersome. 

P(A/ < Xi, I < i < m.Yj < X, Tn < °°) 

= ^ F{Tn = {k+l)/2",C„ = {ki + \)/2'^,Xl<Xi,l<i<m,r <x) 

= ^ P (k/2" < T < {k+\)/2’^,k/r < C\^^X)lin < l)/2", 

^t+(ki+l)/I'’ “2f(^i+l)/2" < 1 < i < '”5h'f+(P+I)/2« -Y{k+l)/2>‘ < 

= Y. P(V2"<T<(^+l)/2”,k,/2”<Ci^+i)/2«<(^/ + l)/2”,l<i<"i) 

X P(X/ < A;, 1 < / < m, Tf < x) 

= P(t„ < oo) P(X/ < Xi, I < i < m,Yt < x). 

As n —)■ oo, the process X' converges to X'■ — X' . We conclude (10). □ 

C^ + ’ ex' 

4. Construction of affine processes on 

In this section we aim at completing the proof of Theorem 1 in the case where the 
process takes values in M™; that is, when n = 0. 

Let X \ ... ,X'”,y be Levy processes satisfying the conditions of Theorem 1 when 
n = 0. Let z G have non-negative coordinates, let Z be the unique solution to (2) 
(when there is uniqueness and let Z be zero otherwise) and let C be the (coordinatewise) 
primitive of Z which starts at zero. (The solution exists thanks to Lemma 3 and is unique 
by Lemmas 4 and 6). 
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For > 0, recall the definition of the multiparameter filtration 

given in (9). Let ^ be the null sets of P and define 


( 12 ) 






= 


f\/ 


Since the processes are independent and the completed filtrations of 

any one of the Levy processes are right-continuous (cf. [Ber96, Prop. 4, Ch. 1]) then 
one can use [CY03, Ex. 2.5] to see that 




‘^'1 


n 




^1 


ti>Si,t>S 


Lemma 7 (Measurability details and the Markov property). 

(1) For any t >0, Q is a multidimensional stopping time: 

.. 

(2) The class 

(13) (f, = {A G ^ : An {c/ < ti,...,cr < tm] e 

is a o-field and the collection > 0) is a filtration satisfying the usual hy¬ 
potheses. 

(3) The following strong Markov property holds: for any t >0, conditionally on 
Q < °°for \ <i <m 




vl vni vm V v\ — fvl V^ 

— Aqi,, — It j — [A ,...,A ) 


and the process on the left-hand side is independent of 

(4) Z is a a > Q)-Markov process. 


Proof. (1) The fact that C is a multidimensional stopping time follows from Lemma 
5 once we note that C is almost surely equal to the maximal solution to (2) 
constructed in Lemma 3 thanks to the fact that solutions have almost surely 
no spontaneous generation (Lemma 6) and the uniqueness of solutions without 
spontaneous generation of Lemma 4. 

(2) It is easy to prove that % is a cr-field. % also contains the null sets JF since 

every contains them by definition. Also, since C if 

s <t, then {fStd > 0) is a filtration. To see that % is right continuous, we only 
need to prove that Ft^fSt = ^s- Let A G Then for t' > t > 5 we have 

An{c/<ti,...,cr<t„}G,F,,„..,,„,,,. 

Since C\ —)■ C( as t ], 5 , we see that 

An{c/ <ti,...,cr<t,4tAn{ci 

as 11 5 and we conclude that 

An{ci<ti,...,cr<t„}G^,,,...,,„„P 

for 5 < t'. Finally, we have already remarked that I t' s, 

proving that A G 

(3) The proof of the Markov type property follows the same pattern as the one in 
Lemma 6. In fact, it is basically the same proof as the strong Markov property 
for Levy processes at a stopping time. 
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(4) First, from Lemma 5 we deduce the existence of a measurable map Ft which 
applied to functions {P,g-^) satisfying H returns the value at t of the (maximal) 
solution h to (8) (which is the unique one when inputting Levy processes plus 
an initial value). Note also that t C/ = — Ci has a cadlag derivative 

t !-)■ Z/ = and satisfies 

zi=Zi + Y,X)j + Y/ 

i ^ 

where 

Yt = Yt+s-Ys. 

Hence Zf+^ = + Since are independent of 

Zl^Ci (which are (f^-measurable), we see that the conditional law of Z^^^ given 
equals the law of Z started at z on the set Z^ — z. □ 

We now consider a martingale which is fundamental to the proof of Theorem 1 . 

Lemma 8 (An exponential martingale). For any u G M'”, the stochastic process M given 
by 

Mt = [F{u) +7?(m) -ZJ ds 

Jo 

is a martingale. 

Proof. Since M has bounded paths on [0,7] for any t > 0, since u G M.'fjt suffices to 
prove that M is a local martingale. Consider the exponential martingale associated to 
any X‘ and to Y: since 

and E(e“'^')=/W, 

and since x i—)■ e“'^ is bounded on E = x M” if w G M™, the stochastic processes 

Ml = - f Ri{u) ds and Nt = - f F{u) e“'^^ ds 

Jo Jo 

are martingales. (Note that the above assertion is true even if Q = oo for some t and i.) 
They are independent since ... ,X'” and Y are independent. 

The random variables [Qd > 0) are stopping times and s <t implies < 

Q. Hence, they constitute a multiparameter time change in the sense of Chapter 6 of 
[EK86]. Consider then the time changed processes 

Mi = M' oCi = _ r zjds. 

Jo 

Problem 19 in [EK86, Ch. 2] tells us that the (multiparameter) time change of the 
m + 1 independent martingales .. .,M"\N gives rise to the m + 1 orthogonal local 
martingales Mf.. .,M",N. Hence 

[M\Mj] = 0= [M\N] 

for all i, j with i j. 

Note that 

j j ' ' 
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Since and are semimartingales whose loeal martingale parts are orthogonal, 

and whose finite variation parts are eontinuous, we ean use integration by parts, the faet 
that eovariation is bilinear and that the eovariation with a eontinuous finite-variation 
proeess is zero (ef. Theorem 26.6.viii in [Kal02]) to obtain 


p“'^'=Loe. Mart. + ^ / Rj{u)Zids + / du 

J ¥j i 


= Loe. Mart.-I- f e“'^¥R{u) ■ + F(u)] ds. 

Jo 


We eonelude that M is a loeal martingale. 


□ 


We deduee the following result, whieh is important in our proof of stability of the 
multiparameter time ehange equation. Indeed, it is important sinee addition is not eon¬ 
tinuous on the spaee of eadlag funetions (with the Skorohod 7i topology), but it is 
eontinuous when the summands do not have eommon diseontinuities, as is diseussed 
for example in Theorem 4.1 of [WhiSO] . 

Corollary 2. Almost surely, for each j E {I,... ,m} the processes o C', 1 <i<m 
and yJ do not jump at the same time. 

Proof. As shown in the proof of Lemma 8, the proeesses e“'^'°^‘I < i < m and 
are semimartingales for 1 < / < m with zero eovariation. Considering a veetor u all 
of whose eoordinates are zero exeept the j-th whieh equals — 1 , we deduee that the 
semimartingales e~^‘'''°^‘, I < i < m and e~^^ have zero eovariation. Sinee e~^^ is of 
finite variation, we see that 

S<1 

Sinee eaeh summand in the right-most side is negative, we eonelude that X‘J oO and 
yj do not jump at the same time. The same argument applies when eonsidering X^J oO 
and A' J o O if i i' sinee at most one is of infinite variation. □ 

As already mentioned in Seetion 2, there exists a unique funetion u) sueh that 
i//( 0 , u) = u and 

dy/(t,u) „ , . 

^ =Roxfr(t,u). 

We also eonsider the funetion 

^{t,u)= / F o \j/[s,u) ds. 

Jo 

In order to prove that the proeess Z whieh solves (2) when n = 0 is a CBI assoeiated to 
the pair of oharaeteristie exponents R and F, it suffiees to see that Z is a Markov proeess 
(whieh is eovered by Lemma 7), and to prove the following lemma: 

Lemma 9. For any z G M™, and u G 


G{s) =EJe 


,\l/{t-s,u)-Zs+(j>(t-s,u) 


Proof. Let 
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for 5 G [0, t]. We will show that G'{s) = 0 for any 5 e (0, t) which implies that 

and hence finishes the proof. 

To this end, write 


G{s + h)-G{s) 


_jg ^^\j/{t—s—h,u)-Zs-f-h+(j){t—s—h,u) _ ^\l/{t—s—h,u)-Zs+0{t—s—h,u) 

_j_ f ^\ff{t-s-h,u)-Zs+^{t-s-h,u) _ ^\lf{t-s,u)-Zs+(j){t-s,u)\ 


Taking expectations in Lemma (8), we see that 

—E _ ^\jf{t—s—h,u)-Zs+(j){t~s—h,u)'^ 


2 rs+h 

TiL 


Ez \^Fo\if{t-s-h,u)+Ro\if{t-s-h,u)- Z,]) dr 


^ E, oxif{t-s,u)+Ro xj/(t -s,u)- Z,]) . 

On the other hand, we can differentiate under the expectation to obtain: 


-E ^^W^^^h,u)-Zs+(j){t-s-h,u) _ ^\l/{l~s,u)-Zs+(j){t-s,u)^ 

^^xif{t-s,uyz,+^{t~s,u) [Fo\]/{t-s,u)-Z,+Ro(l){t-S, w)]) . 


h 

-)■ 


We conclude that G'(s) = 0. 


□ 


Summary and conclusion of the proof of Theorem 1 when n = 0. Existence for solutions 
to (2) are covered by Lemma 3 and is valid more generally. Uniqueness is then covered, 
through the concept of spontaneous generation, in Lemma 6. Lemma 7 then proves that 
the unique solution to (2) is a Markov process and thanks to Lemma 9 we can identify its 
one-dimensional distributions with those of a CBI process associated to the parameters 
of the underlying Levy processes called R and F. □ 


5. Construction of affine processes on M'|! x R" 

Let Z\...,X'”, Y be Levy processes satisfying the conditions of Theorem 1. Let 
R‘ and F be the characteristic exponents of X' and Y as in Equation (7) and let R = 
■ 5 7?"+'”), where we set 7?' = 0 for m -f 1 <i<m + n. With the first m coordinates 
of these processes we solve ( 2 ) to obtain the non-negative processes Z\. ■ •, Z'” analyzed 
in Section 4. We can then (re)define Z by setting Z'”+“' — z + + Yj and note 

that Z solves Equation (2) when /3 = 0. In this case, we can follow the arguments of 
the case n = 0 presented in Section 4 to see that Z is a Markov process and that its 
one-dimensional distributions are characterized by the computation 

valid for u G R” x iW\ where y/ and (j) solve the Riccati equations 

^\j/{t,u) = Ro\j/{t,u) and ^^{t,u) = F o\j/[t,u) 
at at 

with initial conditions yf{0,u) = u and ^ (0, m) =0. This proves Theorem 1 when /3 = 0. 
Affine processes of this type have been dubbed partially additive in [KRSTl 1] since the 
law of z -|- Z under equals Pg+j whenever z has its first m components equal to zero. 
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We now extend the process Z = (Z^,..., just considered to obtain the full proof 
of Theorem 1. To do this, consider the equations 

zf= Zm+j + O c; + 7“+^' + £ cf ^<j<n 

i=l i=l 

where cf'^ = /q zf ds for m+1 < j < m + n (and Z^’-' = for 1 < 7 < m). If we let 
extend the matrix /3 to be (m + n) x (m + n) by adding zeros at coordinates /, j if i < m 
or j < m, the equations become: 

(14) ZP=Z + C<^I5. 

This is a linear stochastic differential equation driven by Z which, of course, admits an 
unique solution. This is for example contained in [Pro04, Ch. 9§V], where the following 
explicit formula is given: 

zf JZ,. 

Jo 

We first construct an exponential martingale, which takes the place of Lemma 8 . 
Lemma 10. For any u G x iM", the stochastic process given by 


Mf = e“'^^ ~ J ^ 


u-Zf 


F{u) + {R{u)+l5u)-Z^ 


ds 


is a martingale. 


Proof. Lemma 8 can be extended to the case /3 = 0, proving that is a martingale. 
Now, note that 

^u-zf ^ ^u-Z,^cfbu^ 

We now apply integration by parts, noting that since C^/3 m is continuous and of finite 


variation, then 




= 0. We then obtain 


.u-zf f eC^h'^[Z,.R(u)+F{u)]e‘^-^^ds+ f e'^'^^Zs^ue'^^ ds. 

Jo Jo Jo 

F{u) + [R{u)+l5u]-zl 


= Loc. Mart. + 


rt p 
„u-Zt 


Adapting the proof of Lemma 9, we see that 




where 1 //^ and (j)^ satisfy the Riccati equations 

^\j/^ (t,u) = Ro \j/l^ (t^u) + P\j/^ {t,u) and ^^^(t,u) = F o\j/^U,u) 
at ot 

with initial conditions y/^ (0, u) = u and (0, u) = 0. 

We now finish the proof that Z^ is an affine process, thereby proving Theorem 1 in 
the remaining case when nf^O. Since we have already determined the one-dimensional 
distributions of Z^, it remains to discuss the Markov property. Note that 

zf+, = zf + 4+, - Z, + /3 zf dr. 
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Therefore, satisfies the same differential equation as but starting at zf and 
driven by Zs+t — Zf. Recall that the first m coordinates of Z^ equal those of Z. Let 
0) be the filtration defined in Section 4 and with respect to which the Markov 
property of Z holds. Since Z is partially additive, then 


the law of Z/_,_.,..., Z^^., 


m+I 


1 

Z^t T ■ ■ 


't+- 


Z“+" given % is 




This shows that the law of zf,_. given equals the law of Z^ under ¥ p which proves 

that Z^ is an affine Markov process whose infinitesimal parameters we have already 
determined. 


6. Stability analysis of the time change transformation: 

APPROXIMATION AND LIMIT THEOREMS 

In this section, we will give a stability analysis related to the stochastic system (2) 
through the deterministic system (8), aiming at a proof of Theorem 2. For the stabil¬ 
ity analysis we need to consider not only the system (2) but a differential inequality 
that turns up naturally. This differential inequality is analyzed in Subsection 6.1. The 
stability analysis is then performed in Subsection 6.2 which enables us to obtain some 
applications to approximations and limit theorems concerning affine processes in Sub¬ 
section 6.3. 

6.1. A differential inequality. Recall the setting of Theorems 1 and 2. If converges 
to C, it might happen that C/ < Q or C^/ > Q. Hence, we can only infer that 

o c; < liminfX'’^’^ o C'’^ < limsupX'’^’^ o C'’^ < o c/. 

^ / 

The following proposition is useful in determining whether the limit X'’-^’^ o C'’^ exists 
for most values of t. 

Proposition 1. Under the setting of Theorem 1, the associated cumulative population 
C = (C^,..., C™) is the unique non-decreasing and continuous process satisfying the 
differential inequalities 

(15) f o C, + yJ ds < Ci - Ci < f o C, + 7/ ds. 

As a preliminary result, let us see that C itself satisfies both sides of the inequality. 

Lemma 11. Almost surely, ift is such that O is constant on an interval to the right oft 
then X' is continuous at Q. Hence, almost surely, for all t > 0.' 

c/-C/= f oC, + Yi ds. 

Jr i 

Proof The first statement of the lemma is obviously true at t = 0. To handle every t >0 
simultaneously, it suffices to prove that for any rational q, if Tfr,, denotes the first zero 

of Z^ after q or zero, depending on if Z^ > 0 or not, then X' is continuous at C{T^q^). 
This is basically a result of quasi-continuity of the Levy processes involved. 

If T is any stopping time with respect to the filtration (^?, t > 0) defined in (13) then 
{Cj, ...,C'f,T) is a stopping time with respect to the multiparameter filtration 


AFFINE PROCESSES ON R™ x R" AND MULTIPARAMETER TIME CHANGES 22 

defined in (12). This follows simply when T takes values in a diserete set : k G N} 
beeause, by definition of we see that 

{4 < ti, • • • T = ak] = « < ti,• • • ,C™ < n{r = aj e 

When r is a general stopping time, we approximate it by the deereasing sequenee of 
stopping times r„ whieh takes the value k/ 2 " if T G [(k — l)/ 2 ",k/ 2 "). 

Let T equal one of the and note that T is the inereasing limit of r„ where r„ 
is the first time after q that Z-' is below l/n or zero depending on if > 1 /n or not. 
We always have T„ <T. If Z^ > 0 then Tn <T for all n. Recall that is a stopping 
time for the filtration o{X‘',,X‘ : s' >0,s <t,i' i),t >0 defined for each i. Since X’' 
is a Levy process with respect to that filtration, by quasi-continuity, we see that X' is 
continuous at Cy. □ 

Proof of Proposition 1. Denote by C any process satisfying the inequality (15). Re¬ 
call that, from Lemma 3, C is obtained as the limit of C\ where solves Equation 2 
driven by processes strictly bigger than and Yf The simple argument presented 
in Lemma 1 implies that C is bounded above by and therefore C < C. We now let 
T = inf {t > 0 : G < Cf}. By continuity, we see that C = C on [0, t] . Let us suppose that 
T < oo to reach a contradiction. If T < oo, there exist I < i, j < m and Eq > 0 such that 
for 0 < £ < Co we have 

/■T+e . . rT+E 

(16) / X'l^oC\ds< / T'^oClds. 

Jx Jt 

When Zij._ > 0 then C' is strictly increasing to the right of T, implying that C' is strictly 
increasing to the right of T and so (16) does not hold. When Z^_ = 0, T cannot belong 
to the interior (or be the beginning) of an interval of constancy of C'. Indeed, Lemma 11 
would then imply that X‘y o Q = X‘’-^ o Q and that O is constant to the right of T which 
would contradict (16). Hence, O increases on any right neighbourhood of T. However, 
recall that C has no spontaneous generation. This implies the existence distinct indices 
z’o, ■ • • ,4 in {1,... ,m} such that Z‘f > 0, 4 = i and oC‘‘-^ is strictly increasing on 
a right neighbourhood of T. Starting with z’o and using the fact that C = C on [0, l], then 
X'/ i,'; oQi-i is strictly increasing on a right neighbourhood of T for every I and hence 
(16) can not hold either when Z^_ = 0 . □ 

6.2. Stability analysis. The following result deals with stability of the multiparameter 
time changes of equation (2) in the deterministic setting of Section 3. We focus on the 
case n = 0 since our arguments can then handle the non-negative coordinates. Hence, we 
will concern ourselves with equation ( 8 ) not only under changes in and g' for z, j = 
1,..., zn, but also with respect to discretization of the transformation itself. Consider the 
following approximation procedure: given a > 0 , called the span, consider the partition 

4 key, k 0 , 1 , 2 ,..., 

and construct a function = (cj^,..., 4 ) by the recursion 

cj( 0)=0 fory^ l,...,zn. 

and for t G [4-i,4): 

m 

+ i^k- 1 ) + (4- 1 )] ^ • 

i=i 


( 17 ) 
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Equivalently, the funetion is the unique solution to the system of equations 

rt ^ . 

fory = l,...,m. 

•^0 

The stability result is stated in terms of the usual Skorohod 7i topology for eadlag fune- 
tions: a sequenee // eonverges to / if each coordinate converges in the usual Skorohod 
Ji topology. This means that for each coordinate //,/ > 1 there exist a sequence of 
homeomorphisms A /,/>! of [ 0 , oo) into itself such that 

fj — f-’ o Xj and A/ — Id converge to 0 uniformly on compact sets. 

We will also use the uniform J\ topology and which is characterized by: a sequence of 
fi converges to / if for \ < j <m there exists a sequence of homeomorphisms A/of 
[ 0 , oo) into itself such that 

// — ° V and A/ — Id converges to 0 uniformly on [0, oo). 


Theorem 3. Let and be cddldg functions which satisfy hypothesis H, and 

suppose that there exists a unique non-decreasing c which satisfies 

ptm _ ft f" .. . 

(18) / ^/!:^oc'(r)+g-^'(r)Jr <c^'(t)-c^'(5) < / +g^'(r)(ir 

i=l i=l 

for s < t, and 7 = 1,..., m. (In particular, c solves ( 8 ) and has a right-hand derivative 
h.) Let T be the explosion time of c defined by 

T = inf {t > 0 : 37 such that c\t) = oo} . 

If for i,j— 1,..., m, and gj —)■ gj in the Skorohod Ji topology, Oi —)■ 0, and 

Cl is any solution to 



£ f'i^ ° c\ ( [s/Oi]Oi)+g\{[s/Oi]Oi) 


i=\ 


ds 


then Cl ^ c pointwise and uniformly on compact sets [0, t). Furthermore, if p'' o c' and 
fJ’’ o c-^ do not jump at the same time for i j and o c' and g do not jump at the same 
time then the right-hand derivatives Z)+c; converge to h 

(1) in the Skorohod Ji topology if x = °°. 

(2) in the uniform Ji topology ifx<°° and we additionally assume that 
for i = 1 ,.. .,min the uniform Ji topology. 


Theorem 2 follows from Theorem 3 thanks to Lemma 1 and Corollary 2. 

In order to prove Theorem 3 we will first prove a series of lemmata. 

Lemma 12. Under the assumptions of Theorem 3, if (ci(t)J > 1) is bounded for some 
t > 0 then Cl ^ c uniformly on [ 0 ,t]. 

Proof Let M be a bound for ci and let K be an upper bound for (fjj / > 1, z = 1,..., m) 
on [0,M] and (g/,Z > 1) on [0,t] (which exists since f^^ —)■ for ij— l,...,m, and 
g)- Lor any 5 G [0, t] we have that 
(19) 

m 

D+c\(s) = oc\([s/oi]oi)+gl([s/oi]oi)]^ < (m+l)K for any z = l,...,zw, 

z=i 
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implying that the family of functions {c^ : / > 1} has uniformly bounded right-hand 
derivatives (on [0,t]) and starting points. (If O/ = 0 we get the same upper bound for 
D^c\ using the equality £>+c/ = Y!iL\ °c\+ gj.) Therefore {cj : / > IJ— 1,..., m} 
is uniformly bounded and equicontinuous on [0,t]- This in turn implies the same for 
the family {q : / > 1} C C([0,r],R'”) on [0,t]- By the Arzela-Ascoli theorem, every 
sequence (cz^,fc > 1) has a further subsequence that converges to a function c (which 
depends on the subsequence). We now prove that c = c, which implies that c/ —>■ c as 
/ _!. oo uniformly on [0, t]. 

Suppose that 4 is such that has a limit c as k uniformly on [0, t]. Since each 
f'j has no negative jumps, we get 

liminf/'’-^ =/^(y) and limsup/'’-^(jc) =/'’“'(y) 

X^y 


SO that 


oc‘ < liminfo c] and limsup oc\ < o c'. 

fc—k k k k 


Using Patou’s lemma we get 


rt ^ 


\^f-°c\r) +g^_{r)]^dr < c^{t) —c\s) 


/=1 


< / °c\r)+g^{r)]^dr 

Js i=i 


for each j = 1,... ,m. 


□ 


Proof of Theorem 3. For the convergence of the ci to c, we will argue along sequences 
4 —^ considering the following two cases: (c 4 (t)) is bounded or goes to oo. The 
former alternative is handled by the previous theorem, for the latter we prove that —)■ c 
pointwise on [0, t]ask^ The conclusion is that ci^c pointwise on [0, oo) and hence, 
by the previous lemma, uniformly on compact sets of [0, t) . 

Suppose that 4 —^ is such that ||cz^ || ^ —)■ oo. For any .r > 0, consider the sequence 
C 4 A T = [c]^ A T,..., A jc). Note that it is uniformly bounded. Let ^ be a common 
bound for on [0, x] and for g^ on [0, t]. For any ^ G [0, t] 


D, 




1=1 




< (m+l)i^, 


so that the sequence Av is uniformly bounded and equicontinuous on [0,t]. Let c be 
its uniform limit on [0,t]. If c^{s) < x for all j = 1,... ,m we can argue as in the proof 
of the previous lemma to see that c = c on [0, 5 ]. If c^s) > X for any j = 1,..., m, we 
see that both d and both reach x at the same point s' <s and hence c^{s) = c^{s) A v. 
Hence A v —)■ d A v. Since v is arbitrary, we see that —)■ d pointwise on [0, t ], even 
if |c(t)| = 00 . 

Let hi = D+c/ = (L>+c|,..., L>+c™) and h = D+c = (L>+c^,..., D+c™). We now prove 
that hi ^ h in the Skorohod Ji topology if the explosion time T is infinite. Recall that 
h = Y 4 L 1 f‘ oc‘ +g and that when C/ = 0 then hi = YpLi fl °c‘i + gi while if a/ > 0 then 
kjit) = [YJiLifi’'’°c\{[s/ai\ai) + gj{[s/oi\oi)]^. Assume that a/ = 0 for all I (the 
arguments are analogous when < 7 / > 0 ), then the assertion hi ^ his reduced in proving 
that : fl oc‘ —)■ /' o c' for all / = 1 ,..., m, which is related to the composition mapping 
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on Skorohod space, and then dedueing that Y 4 L 1 fl °c] +81 ^ E^i /' o c' + g, whieh is 
related to eontinuity of addition on Skorohod spaee. Both eontinuity assertions require 
eonditions to hold: the eonvergenee // oc\ —>■ /' o c' ean be dedueed from [Wu08] if 
we prove that /' is eontinuous at every point at whieh is diseontinuous, and the 

eonvergenee // o c| + g/ —)■ T.T=\ +8 will hold beeause of [WhiSO, Thm. 1.4] 

sinee we assume that /' o c' and fj o will not jump at the same time nor as the same 
time as g. Henee the eonvergenee hi ^ his redueed to proving that /' is eontinuous at 
the diseontinuities of If c' is strietly inereasing then is eontinuous. When 

c is not strietly inereasing, we will use the assumed uniqueness of (18) to prove that /' 
is eontinuous at the diseontinuities of The proof eonsists in two steps. First we 

will prove that /" is eontinuous at c‘{s) and the we will use this faet to prove that the 
rest of the eomponents of f is eontinuous at the same point. 

We know prove that is eontinuous at the diseontinuities of Suppose that 

is diseontinuous at t. Let ^ ) and t = (c')~^(t). Then c' is eonstant 

on on [5,1] while c' < Ton [0,5) andc' > von Sinee D+c' = Y!f=if +g‘ = 0 
on [ 5 , 1 ), we see that + 8 ' is eonstant on [ 5 , 1 ). We assert that 

m 

inf{y>0:/’'(y) = - ^ f'’^oc\s)+g\s)]=x. 


Indeed, if reaehed —E™=i +g'(^) at v' < x, there would exist s' < s 

sueh that 


f-'oc\s)+ £ /•‘oc'\s)+g\s)^ 0 >roc\s')+ £ 


SO that aetually we have that — E/Li +g'(5') is eonstant in [s',t). Henee, c' 

has spontaneous generation, whieh implies that there are at least two solutions to (18): 
one that is eonstant on ( 5 ', 5 ), and cf This eontradiets the assumed uniqueness to (18). 

Having proved the eontinuity of at x, we need to prove that for eaeh j 7 ^ i, that /'’-i 
is eontinuous at v. To this end let us reeall that c\s) = x, and eonsider a sequenee {x„} 
sueh that v„ t v, therefore there exists another sequenee { 5 ,,} sueh that Sn t So using 
the eontinuity of at v we have that 


0 < lim (/'’'oc'(5„) + y /^’'oc-''(5„)+g^'(5„) j 

m 

<f'‘oc‘{s)+ Y, f'‘ocJ{s)+gJ{s)=0. 

Therefore we obtain that 

m m 

lim Y f''°<^'’i^n)+8^isn)= Y f’'°<^'’i^)+8''{s). 


And henee we ean eonelude, using that the funetions /i’' and gj, are non-deereasing 
that 

lim fj'' o c^{sn) = o c-i(j') 

n—^00 


for all j 


□ 
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6.3. Applications of the stability analysis. In this subsection, we apply the stability 
analysis of Subsection 6.2 to give a proof of Corollary 1. 


Proof of Corollary 1. We first analyze the action of scaling on Since 


- 


rait 


d m 

^ + Lyr,j,ic 




'W 


ds 



ibj^ bjh V 


ds, 


we see that {clft/bj,t > 0,1 < y < m) is the Euler type approximation of span 1/a/ 
applied to > 0,1 <i,j<m). Note that by hypothesis, the span 1/a/ 

goes to 0 as / —)■ oo. Also, the right-hand derivative of Cflt/b\ equals Z/’^az/Z?/. 

Also, we have assumed {X^'^'\b\t)],b\/ai,t > 0,1 < y < m) converges to A'’’. Since 
b\/a\ —)■ 0, we see that A'’' is spectrally positive. If i f j then A'’-' is a subordinator. 
(W^e have only assumed convergence in the Skorohod 7i topology. However, the same 
arguments as in the proof of Corollary 7 of [CPGUB13] gives us convergence in the 
uniform topology in case of explosion.) Finally, since A^’’,... ,A™’' are independent, 
we are in position to construct the CB process Z starting at and constructed from A 
and y = 0 in Theorem 1. From Theorem 2, we deduce that (Z'/^fai/bj) converges to Z 
as Z —)■ oo, which proves Corollary 1 . □ 
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